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Abstract 

The classical stochastic model of cosmology recently developed by us is 
reconsidered. In that approach the parameter w defined by the equation 
of state p = wp was taken to be fluctuating with mean zero and we com- 
pared the theoretical probability distribution function(PDF) for the Hubble 
parameter with observational data corresponding to a universe with matter 
and vacuum energy. Eventhough qualitative agreement between the two was 
obtained, an attempt is herein made to introduce a more realistic assump- 
tion for the mean of w and use it for the calculations. In the present theory 
the mean values of both p and w are taken to be nonzero. The theoretical 
and observational PDFs are compared for different epochs and values of the 
Hubble parameter. The corresponding values of the diffusion constant D 
obtained are approximately constant. We use the scatter in the observed 
redshift-magnitude data of Type la supernova to place limits on the stochas- 
tic variation in expansion rate and consequently, on the stochastic variation 
of the equation of state. 
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1 Introduction 



The uncertainty in the determination of the Hubble parameter H, which is a 
measure of the expansion rate of the universe, is one of the most intriguing 
issues in the history of cosmology. The origin of the uncertainty is obvious 
from the redshift-magnitude diagram (Hubble diagram); despite rigorous at- 
tempts to control the random errors in measurement, there is a clear scatter 
in it, though it is now possible to narrow down this to a great extent. But, 
as we show later in this paper by using the redshift-apparent magnitude data 
for Type la supernovae (SNe la) [1], the scatter increases as we go to higher 
redshifts. In a recent work [2], we have attempted to explain this scatter as 
arising from an inherent stochastic or nondeterministic nature of the Hubble 
parameter. It was shown that a fluctuating w-factor in the equation of state 
p = wp will lead to this kind of behaviour for H. In order to arrive at the 
notion of a stochastic equation of state, we consider the following. Several 
recent measurements [1, 3, 4, 5, 6] indicate that the present universe contains 
components other than ordinary matter and radiation, like vacuum energy, 
quintessence etc.. Let the equation of state for each component, with den- 
sity pi, be written as Pi = WiPi, (i = 1, 2, ..), where Pi is the corresponding 
pressure [yji = for dust, 1/3 for radiation, —1 for vacuum energy, etc.; in 
general, — 1 < Wi < +1). The total energy density is p = J2iPi an d the total 
pressure p = J2iPi- From the law of conservation of total energy-momentum 
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tensor expressed in the form 



3-p(l + w) 



a 



(1) 



which follows from the Einstein equation and one can write the total w-factor 

as 



may be expressed as p\ = —3(a/a)pi(l + Wi), is only an extra assumption 
since it does not follow from the Einstein's field equation. Equivalently, it can 
be stated that in a many component fluid as in the above case, the Einstein 
equations, along with the equations of state of individual components are 
insufficient to determine the individual pi's. Thus it is more general not to 
assume conservation of individual components and this will lead to creation 
of one component at the expense of other components. Since they are not 
uniquely determined by the field equations, such creation can be considered 
as sporadic events, like those occurring in galactic nuclei, which can result 
in fluctuations in the ratio pi/ p. Here we make the assumption that this will 
lead to a stochastic equation of state, as seen from equation (2) above. Con- 
sequently, also the expansion rate will be fluctuating and the equation for 
the Hubble parameter will appear as a Langevin type [7] equation. Recently 
Padmanabhan [8] has proposed a similar scenario in which the observed cos- 




(2) 



where a is the scale factor. The conservation of individual components, which 
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mological constant arises from quantum fluctuations of some energy density; 
consequently the FRW equations are to be solved with a stochastic term on 
the right hand side. In [2], it was assumed for the sake of simplicity that 
space sections are flat and w is a Gaussian 5-correlated stochastic force with 
zero mean. With these assumptions, we have written the Fokker-Planck (FP) 
equation, whose solution gives the theoretical probability distribution func- 
tion (PDF) for H Q at time to, denoted as W(h,t ) (where H = lOO/i km 
s _1 Mpc -1 ; the subscript denotes the present epoch). Using the redshift- 
apparent magnitude data /i for SN la used in [1], we computed the obser- 
vational PDF p(h\fi ) for h in the present universe, again assuming its space 
sections to be flat. This PDF arises from the point-to-point variance of the 
Hubble flow. We compared the two plots for the present universe and found 
them to agree well, for a value of the diffusion constant D appearing in the 
FP equation equal to 4 x 10 13 s. 

This result is a first step towards an understanding of the anomalous 
scatter in the Hubble diagram at high redshifts. However, there are certain 
refinements to be made in our analysis. One drawback of the above scheme of 
comparing these two PDF's is that when we derived W{h, t), the assumption 
was made that w has mean value zero, whereas the observational PDF p(h\fj, ) 
was evaluated for a model which contains dust and vacuum energy, which has 
the mean total pressure negative. Instead, if we had used in this evaluation 
the expression for the distance modulus for a flat universe which is dust 
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dominated (i.e., with w = 0), an observational PDF will be obtained, but 
the best fit value for h would be ridiculously low. (This is all the fuss about 
the new observations - that they are incompatible with an JIa = flat model.) 

Another shortcoming is that though in both cases we take the PDF for 
h, it remains to be explained how legitimate is the comparison of W(h,t ) 
for the present universe with a PDF p(h\fi ) evaluated using the data that 
include high redshift objects, which belong to the distant past. 

In this paper, we attempt to rectify these two defects and to make a 
more rigorous test of the stochastic assumptions using observational data by 
(1) comparing both the theoretical and observational PDF's evaluated for 
the same model, which is an alternative flat model [9] and (2) evaluating the 
observational PDF p(hj \fi oj ) for the Hubble parameter at the same epoch tj as 
that in the theoretical PDF W(hj,tj). This procedure helps us to compare 
the theoretical and observational PDF's for the Hubble parameter for the 
same model, and at the same epoch. The value of the diffusion constant 
evaluated at any time is obtained as nearly a constant, in agreement with 
our assumptions. A novel feature in our present approach is that we evaluate 
the observational PDF for Hubble parameter at various instants in the past, 
also with an objective of justifying our assertion that the scatter increases as 
we go into the past. 

The paper is organised as follows. Section 2 gives a brief review of the 
alternative model to be used and then develops the stochastic approach in 
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it. Section 3 gives the new technique of finding the PDF for h at any time in 
the past or present epochs and the results obtained on comparison between 
the theoretical and observational PDF's. The conclusions derived from it are 
discussed in section 4. 

2 Stochastic approach in the new model 

In all FRW models, the Einstein equations, when combined with the con- 
servation of total energy-momentum tensor can be written in terms of the 
Hubble parameter as 

H = -H*-^(p + 3p). (3) 

Overdot denotes time derivative. If we restrict ourselves to flat models, then 
(with p = wp), 

H = -\h\1 + w). (4) 

In [2], we considered this flat case and assumed that w is a Gaussian 
5-correlated Langevin force term with mean value zero. This means that the 
mean total pressure of the universe is zero, the same as that for dust. But 
many recent observations are incompatible with this model and hence, as 
mentioned in the introduction, we look for a more observationally correct, 
but simple model to apply our stochastic treatment. 

The deterministic model [9] we propose to use is the one in which the total 
energy density obeys the condition p + 3p = and hence having a coasting 
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(a oc t) evolution. This model is derived on the basis of some dimensional 
considerations in line with quantum cosmology. If we assume that the energy 
components in this model are dust and vacuum, then the above condition 
gives Pm/p v = 2 and if they are only radiation and vacuum, then p r /p v = 
1. In [9], it was shown that in this model, most outstanding cosmo logical 
problems like those of flatness, horizon, monopole, entropy, size and age of 
the universe, and the cosmological constant are absent. It was also shown 
that this model can solve the problem of generation of density perturbations 
at scales well above the present Hubble radius and that it can generate such 
density perturbations even after the era of nucleosynthesis. Though it is 
mentioned in the paper that recent observations on the redshift-apparent 
magnitude relation are at variance with the predictions of this model, in a 
recent communication [10], this issue was studied in detail and found that 
those observations do not provide any strong evidence against it. 

It may be noted that the underlying model described above deviates from 
'main stream' cosmology to some extent and the most significant cosmological 
observations such as nucleosynthesis, large scale structures and microwave 
background radiation, in this model are not well-studied. However, to repeat 
our statements above, we choose this model to implement our stochastic 
approach mainly due to its simplicity and its capability to explain the recent 
supernova data. 

In view of the fact that this model has w = —1/3 in the deterministic 
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case, we rewrite the above equation (4) with w' = w + (1/3), as 



2 V3 

We now assume that w' fluctuates about its zero mean value and is 5- 
correlated. Making a substitution 

1 



the above equation becomes 



X =H 



3 

x = 1 + - w '. (6) 



When w' = 0, this equation is analogous to that of a particle moving in a 
medium with constant velocity. With a fluctuating wf, the analogous particle 
is subjected to random forces as it moves. We write the FP equation for the 
distribution function for the variable x as [7] 



dW' . 



- D m^L + D m ai 



dx dx 2 



W'(x,t), (7) 



the solution of which gives the PDF W' for the variable x at time t. Here 
the drift coefficient is the constant velocity term, equal to unity, and the 
diffusion coefficient = D is assumed to have some constant value, to be 
determined from observation. The FP equation is solved by first assuming 
that the variables are separable; 

W'(x,t) =<p n (x)e- x " t , (8) 



where (f> n (x) and \ n are the eigen functions and eigen values of the Fokker- 
Planck operator 



(x) H 7T^—{X) 



(9) 



dx dx 2 

Now we define two more functions in order to get a solution for the FP 
equation; 



$ (x 



X 



DW dX ~ D 



and 



ip n (x) = exp <j) n (x) = exp ^- 



[X 



(10) 



(11) 



2/ r v ' L \ 2D, 
where $ (x) is treated as a stochastic potential and ip n (x) are the eigen 
functions of the Hermitian operator L defined as 



L = exp ^— j L FP exp (- 



(12) 



Making use of equations (8) and (11), the time independent part of FP 
equation becomes 



d 2 A 
dx 2 



(x) 



1 A r 



AD 2 D 



i>n (x) = ~k 2 i) n (x) 



Here, 



k = ± 

The most general solution to (7) is 



D 4D 2 



(13) 



(14) 



W (x,t) = ^c n e- A "Vn (x). 

n=0 



(15) 
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For \ n < jp, the solution ip n (x) is exponentially diverging, which is not an 
admissible solution. Thus we conclude that A n > ^ so that 

(f) k (x) = A exp (Jjj + ikx^j . 

We make use of the completeness relation 

/+oo 
^* k {x)^ k {x')dk 
-co 



(16) 



(17) 



to specify the initial condition x — x' at t — t'. The transition probability 
for the variable to change from x' at time t' to x at time t is [7] 



P(x,t\x',t') = e LFP ^5(x-x') 



(18) 



cxp 



/+oo 
dk^l{x)i) k {x') exp [-X(k)(t - t')} 
-oo 



2JnD(t - t') 



exp 



[x - x) -(t- t') 



i 2 



4D(t - t') 



where we have chosen A = for normalisation purpose. For the special 
initial condition W'(x, t') = 5(x — x'), the transition probability P(x, t \ x', t') 
is the distribution function W'(x,t) [7]. In our case, we have the initial 
condition x = x' = at t = t' = 0, so that 



W (x,t) = P(x,t\0,0) 



1 



■ exp 



{x-ty 

ADt 



(19) 



which is in Gaussian form and has its peak moving along in such a way that 
the expectation value of the variable is (x) = t. This corresponds to the 
deterministic solution of (6). The width of the Gaussian is found from the 
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va.na.neo a~ = ((x— (x)) 2 ^ = 2Dt and we find a > (x) till t = 2D. We 
can immediately rewrite this distribution function in terms of the stochastic 
Hubble parameter H as 



W"(H,t) 



1 



: CXp 



my 



AH 2 Dt 



(20) 



2H 2 ^D~t 

With H = 100 h km s" 1 Mpc" 1 , t = t 17 x 10 17 s and D = D 17 x 10 17 s, the 
PDF W(h,t) can be written as 



W(h,t) 



3.0856 



: exp 



(3.0856 - htvr) 
WD 17 t 17 



2l 



(21) 



2h 2 ^/nD 17 t 17 

For the range of values of interest, 1 < t 17 < 5 and Dn ~ 10 -3 , W(h,t) is 
approximately a Gaussian. For fixed D, the half width of the Gaussian is 
found to increase as we go to lower values of t. 

3 PDF for H from observational data 



From the Hubble diagram, one can find the PDF for the present Hubble 
constant H in the following way. The traditional measure of distance to a 
SN is its observed distance modulus \x Q = rriboi — Mm, the difference between 
its bolometric apparent magnitude and absolute magnitude. In the FRW 
cosmology, the distance modulus is predicted from the source's redshift, z, 
according to 



H P = 5 log 



D, 



+ 25, 



IMpc 

where D L is the luminosity distance, found as 

D L =r ja (t )(l + z). 



(22) 



(23) 
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r° cat 

r ' = L m- (24) 



a(to) is the present scale factor and rj is the radial coordinate of the SN la 
which emitted the light at some time tj in the past. In flat FRW models, r,- 
is found as 

r*o c dt 

For the coasting model discussed in the previous section, for k — 0, rj can 
be evaluated as 

r,- = -^ r* = -^ln(l + *), (25) 
3 a(t )J tj t a(t ) K h K J 

so that 

^ = ^±^(1 + 4 (26) 

-no 

One can use this expression in (22) to obtain the predicted distance modulus 
of an object with redshift z. Conventionally, assuming that the observed 
and predicted distance moduli coincide, one can find a value of H . For 
a collection of objects, also one can find the likelihood for H , from a x 2 
statistic; 

x 2^ (/^-/V) 2 , (27) 
i a i 

where <7j is the total uncertainty in the corrected peak magnitude of SN la. 
For the special model we are considering, h is the only parameter and the 
normalised PDF can now be obtained as [3] 

(h\ \ ex P(-X 2 /2) , . 

P < fc l"»» = /»*e X p- x V2' (28) 

As in [2], we compute p(h\/i ) for the new model using the SNe data in 
[1], which corresponds to their Fit C and attempt to compare p(h\/i ) with 
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the PDF W(h,to) to evaluate the diffusion constant D appearing in this 
expression. It is found that the two curves, shown in figure 1, coincide for 
a value of D m 2.36 x 10 13 s. (This corresponds to an age 4.8583 x 10 17 s.) 
Since our primary objective is to make an order of magnitude evaluation of 
D, we chose a fiducial absolute magnitude for SN la in computing \i Q , equal 
to —19.3 mag. Slight variations in this quantity will not significantly affect 
D, though the best fit value for h may change. 

In the above, we compared the theoretical and observational PDF's for 
the same alternative model and thus it does not have the first shortcoming 
mentioned in the introduction. The other incompatibility which still exists 
can be explicitly stated as follows: W(h,t ) is the PDF for the Hubble pa- 
rameter of the present universe and it contains the diffusion constant D . 
But p(h\fi ), which we attempt to identify with W(h, to), depends on the 
scatter in the Hubble diagram for all ranges of redshift. For instance, if we 
include more high redshift objects in our sample, the scatter would be larger 
and hence the half-width of the distribution p(h\/i ) will be larger. This, in 
turn, will affect the computed value of D, which is quite unreasonable. 

This problem can, however, be overcome if we agree to compute p(hj\/i j) 
for each value of redshift z (or for small enough redshift intervals centred 
about such values), and compare these with W(hj,tj) that corresponds to 
the same epoch tj. To do this, we modify (26) by re-evaluating r,- in (25) in 
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a different way. One can also write, for the new deterministic model 

ctj [to dt c , . , . 

fj = — I — = — — -lnl + A 29 

J a(tj) J tj t Hja(tj) y ' y ' 

so that 

D L = C{X t zf ln(l + s). (30) 

Evaluating \i p using this expression in (22), we can calculate y 2 and hence 
also p(hj\/i j), which is the PDF for the Hubble parameter at some particular 
value of z. We divide the data in [1] for various ranges around z = 0.05, 0.15, 
0.35, 0.45, 0.55 and 0.65, each with Az = 0.05. The PDF for the average 
Hubble parameter for such intervals is calculated with an expression identical 
to ( 28). The results are plotted in figure 2 along with the corresponding 
W(hj,tj) which overlaps with them. The relevant parameters are given in 
table 1. 

4 Conclusions 

It can be noted from figure (2) and table 1 that for the intervals with larger 
values of z, the 68.3% credible region of p(h\fj, j) has a half-width ah which 
also increases (The intervals with centre at z — 0.15 and 0.35 are exceptions 
to it, but this may be due to the fact that these intervals contain only very few 
objects. As more SN la are observed in these redshift bins, an accurate pic- 
ture will emerge.). This behaviour is the one expected from theory, as noted 
while plotting the theoretical PDF (21). The value of the diffusion constant 
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D evaluated for various intervals, however, does not show any dependence 
on z. This justifies our assumption (taken for the sake of simplicity) that D 
is some constant. 

As remarked in the introduction, a novel feature of the present analysis 
is that we compute the PDF for H at various epochs in the past. However, 
there is a limitation, even in the present analysis; the intervals we consider 
are with Az = 0.05 and this value may not be small enough to give correct 
answers. The procedure can be improved very much in the future, when the 
number of observed SN la becomes large. 

It should be kept in mind that the measurements are made with finite 
accuracy and hence have a scatter arising from intinsic statistics of error. 
Here we have not shown that the scatter in the Hubble diagram is signifi- 
cantly more than that expected from the known errors. The present theory 
can survive only if p(h \ fj, ) does not become narrower with increase in ac- 
curacy of measurement. Until it is confirmed by future observations like the 
Supernova Acceleration Probe (SNAP) [11], which aims to give precise lu- 
minosity distance of 2000 SN la upto z — 1.7 every year, the computed 
value of the diffusion constant can not be claimed as accurate. Given this 
circumstance, it is safer to conclude for the present that the scatter in the 
SN la data places limits on the stochastic variation in the expansion rate and 
consequently, on the stochastic variation of the equation of state. In other 
words, it is justifiable to consider the value of D computed by us using SN 
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Ia data as an upper limit to a possible diffusion constant for the universe. 
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Table 1 



Redshift 


No. of SNe 


Best fit 


Standard 


Age in 


Diffusion 


z 


in the 


value of 


deviation 


units of 10 17 s 


constant 




interval z ± 0.05 


h 


Oh 


in 


D s 


0.05 


15 


0.693 


0.011 


4.4502 


0.5775xl0 14 


0.15 


3 


0.772 


0.025 


3.9987 


2.147 xlO 14 


0.35 


5 


0.830 


0.024 


3.7204 


1.5 xlO 14 


0.45 


15 


0.875 


0.017 


3.528 


1.655 xlO 14 


0.55 


7 


0.985 


0.026 


3.1345 


1.12 xlO 14 


0.65 


6 


1.043 


0.030 


2.959 


1.226 xlO 14 
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Figure caption: 

Figure 1: Observational and theoretical PDF's vs h, using the apparent 
magnitude- redshift data for Type la supernovae as given in [1], which cor- 
responds to their Fit C. The continuous line is for the observational PDF 
whereas the dotted line gives the theoretical PDF. 

Figure 2: Observational and theoretical PDF's vs hj (where hj corresponds to 

the Hubble parameter for the epoch centered about redshifts z = 0.05, 0.15, 0.35, 0.45, 0.55, 0.65), 

using the apparent magnitude-redshift data for Type la supernovae as given 

in [1], which corresponds to their Fit C and which lies in the interval z ±0.05. 
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